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Abstract 

We consider a diffusive Coupled Map Lattice (CML) for wliicli tlie local map 
is piece-wise affine and has two stable fixed points. By introducing a spatio- 
temporal coding, we prove the one-to-one correspondence between the set of 
global orbits and the set of admissible codes. This relationship is applied to the 
study of the (uniform) fronts' dynamics. It is shown that, for any given velocity 
in [—1,1], there is a parameter set for which the fronts with that velocity ex- 
ist and their shape is unique. The dependence of the velocity of the fronts on 
the local map's discontinuity is proved to be a Devil's staircase. Moreover, the 
linear stability of the global orbits which do not reach the discontinuity follows 
directly from our simple map. For the fronts, this statement is improved and as 
a consequence, the velocity of all the propagating interfaces is computed for any 
parameter. The fronts are shown to be also nonlinearly stable under some restric- 
tions on the parameters. Actually, these restrictions follow from the co-existence 
of uniform fronts and non-uniformly travelling fronts for strong coupling. Finally, 
these results are extended to some C°° local maps. 
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1 Introduction 



The occurrence of interfaces is a very common phenomenon in extended systems 
out of equihbrium. These patterns appear, in systems with many equihbrium 
states, as solutions hnking two different equihbria, namely two phases When 
the time evolution of the system shows the motion of an interface resulting in 
an invasion of one phase into the other, the phenomenon under consideration is 
usually called a front. These structures arise in various experimental systems 
such as reaction-diffusion chemical systems 0, alloy solidification ||^ and crystal 
growth Q. 



From a mathematical point of view, the description of the front dynamics 
using space-time continuous models is now complete, at least for several one- 
dimensional Partial Differential Equations [^]. Also, in discrete space and contin- 
uous time systems, i.e. in Ordinary Differential Equations, the dynamics of these 
interfaces is well understood 0. However, in space-time discrete models such 
as the Coupled Map Lattice (CML), the study is not as complete]^; except for the 
one-way coupled model 0. CML's have been proposed as the simplest models 
of space-time discrete dynamics with continuous states and they serve now as a 
paradigm in the framework of nonlinear extended dynamical systems [10[ . 



The phenomenology of the interface dynamics differs between discrete space 
and continuous space systems. In the former case, varying the coupling strength 
induces a bifurcation^ from a regime of standing interfaces to a regime of fronts in 
which the propagation can be interpreted as a mode- locking phenomenon |^. 
These effects are assigned to the discreteness of space and are well-known in 
solid-state physics ITT . 



In previous studies of the piece- wise affine bistable CML, the problem of 
steady interfaces was solved and the fronts' structure was thoroughly investigated 
by means of generalized transfer matrices [|12| , |13|. Nevertheless, this technique 



did not allowed us to prove the existence of fronts of any velocity, nor to under- 
stand clearly their dependence on the parameters, in particular that of the fronts' 
velocity. 



In the same CML, using techniques employed for the piece-wise linear map- 



pings IIJ, we now prove the one-to-one correspondence between the set of 
orbits that are defined for all times, i.e. the global orbits, and a set of spatio- 
temporal codes (Section 2 and 3). Our model is locally a contraction, hence 

■^The problem of travelling waves has been investigated in other space-time discrete models 
such as the chain of diffusively coupled maps for which the coupling is different from that in 
CML's II . 

*This bifurcation is of saddle type and is accompanied by a symmetry breaking p2|. 
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these orbits are shown to be hnearly stable when they never reach the local 
map's discontinuity (Section 4). Further, using the orbit-code correspondence, 
the existence of fronts is stated and their velocity is computed (Section 5). In 
Section 6, the linear stability of fronts with rational velocity is given for a large 
class of initial conditions. In the following, we study the dynamics of the propa- 
gating interfaces. In particular, their velocity is computed for all the parameters 
(Section 7). Using these results, the nonlinear stability, i.e. the stability of fronts 
with respect to any kink initial condition, is proved in Section 8 using a method 
similar to the Comparison Theorem P]. This result holds for any rational ve- 
locity provided the coupling is small enough. We justify such a restriction by 
the existence of non-uniform fronts for large coupling which co-exist with the 
fronts. Finally, some concluding remarks are made, in particular we emphasize 
the extension of these results to certain C°° local maps. 



2 The CML and the associated coding 

Let M e M be fixed and M = [-M, M]^ be the phase space of the CML under 
consideration. The CML is the one-parameter family of maps 

: M — >M 

where = {x\} is the state of the system at time t. This model is required 
to be representative of the reaction-diffusion dynamics. Therefore the new state 
at time t -|- 1 is given by [ IC ] : 

= (F,x*). = (1 - e)/(x*) + I {f{xU) + /(4+i)) e Z. (1) 
Here the coupling strength e G (0, 1] and we choose the local map / to be bistable 



and piece- wise affine |jT3|: 



ax + {1 — a)Xi if x < c 
ax + {1 — a)X2 if x ^ c. 



where a G (0, 1) and —M ^ Xi < c ^ X2 ^ M. These conditions ensure the 
existence of the two stable fixed points Xi and X2, the only attractors for /. 

This local map reproduces qualitatively the autocatalytic reaction in chemical 
systems 0, or the local force applied to the system's phase in the process of 
solidification @. 

For the sake of simplicity we assume Xi = and X2 = 1. This is always 
possible by a linear transformation of the variable. 
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For a state x*, the sequence 6** = {91}^^^ defined by 



sG 

9i ^ 



if xl < c 

1 if x\-^ c, 



is called the spatio-temporal code or the code unless it is ambiguous. 



3 The orbit-code correspondence 

The study of the orbits, in particular those that exist for all the times t E Z, can be 
achieved using their code. In this section, we first compute explicitly the positive 
orbits^ for any initial condition. Then we prove the one-to-one correspondence 
between the global orbits and their code. 

Notice that the local map can be expressed in terms of the code 

fixl) = ax\ + (1 - a)ei. 

By introducing this expression into the CML dynamics, one obtains a linear non- 
homogeneous finite-difference equation for x\ in which the code only appears in 
the non-homogeneous term. Using the Green functions' method, this equation 

may be solved and the solution, as a function of the code \ 9^\ and of 

L ^ J sGZ,to^fcs;t-i 

the initial condition is given by 



t—to k ^ t—to 

^l=Yl 5Z ln,kOlZn + J2 ^n,t-toxl% Vt > to and s e Z, (2) 

— — 1 — (I — 

k=l n=—k n=— (<—<()) 



where the coefficients ln,k satisfy the recursive relations 



and 



ln,k = (1 - a)Si^k { (1 - e)Snfi + ^{Sn,i + Sn-i) ) Vn G Z and /c < 1,[ 



ln,k+i - a(l - e)ln,k - Y iln+i,k + ln~i,k) = Vn G Z and k ^ 1. 



From the latter, it follows that for e G (0, 1) 

ln,k = iff \n\ > k OT n = k = 0, 
and, one can derive the bounds 

^ ln,k ^ {I - a)a''-\ (3) 

^i.e. the orbits for i ^ 

^5n,k is the Kroenecker's symbol 
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and the normalization condition 

E ln,k = 1. 

n,fcez 

Further properties of these coefficients are given in Appendix ^ 

The study is now restricted to the global orbits {x^}t£z- That is to say, we 
consider the sequences which belong to 

A = [{x'}tez e : F,(x*) = a;*+^ Vt G z} . 

For such orbits, taking the limit ^ —oo and using the bounds (H) in @) leads 
to the relation: 

+00 k 

^l = T.i: ^nA'-n yS,tEZ, (4) 
k=l n=~k 

which gives the correspondence between these orbits and their code as we state 
now. 

By the relation (^, all the orbits in A stay in [0, 1]^. Hence their code can be 
uniquely computed by 

Ol=U-c\+l ys,tez, (5) 

where [.J stands for the floor function^. 

From (^) and it follows that for any orbit in A, its code must obey the 
following relation: 

+00 k 
k=l n=—k 

which is called the admissibility condition. Then, conversely to the preceding 
statement, for a sequence 9 G {0, 1}^ that satisfies this condition, there is a 
unique orbit in A, given by (Q). 

The spatio-temporal coding is an effective tool for the description of all the 
global orbits of the piece-wise affine bistable CML [|I6|. These orbits are impor- 
tant since they collect the physical phenomenology of our model. 

i.e. [xj e Z and x — 1 < [xj ^ a; 
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4 The linear stability of the global orbits 



We prove in this section the stabihty of the orbits in A with respect to small 
initial perturbations. Let the norm 



||X|| = sup \Xs\, 

for x G and 

(Lx)^ = a (^(1 - e)xs + ^{x^+i + x^-i) j Vs G Z, 

be the CML linear component. Notice that L is invertible on if e < | and 

IlLxll 

L = sup -r. — TT- = a. 

XytO \\X\\ 

The linear stability is claimed in 

Proposition 4.1 Let {x^j^^j^ be an orbit such that 

36 > : |x* - c| > 6a* Vt ^ 0, s G Z. 
Then for any initial condition in a neighborhood of , i.e. 

we have 

lly* < (5a* Vt ^ 0. 

Equivalently, x* and have the same code for all times. 

Proof: The relation (|1|) with the present local map can be written in terms of the 
operator L 



= Lx* + ^L^* (6) 
a 



where 6** = {9\}s&i. is the code of x*. Using this relation, one shows by induction 
that the codes of the two orbits remains equal for all the times; also using @, 
the latter implies the statement. □ 

Notice that this assertion is effective for the orbits in A that never reach c, 



since in this situation 6 can be computed (see Proposition |6.1| below for the case 
of fronts). Further, because our system is deterministic, when this proposition 
holds for an orbit {x*} and a given initial condition y, it cannot hold for a different 
orbit {x*} and the same initial condition unless both these orbits converge to 
each other. Hence, using this statement, one may be able to determine the (local) 
basin of attraction for any orbit in A that never reaches c. 
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5 The existence of fronts 



We now apply the orbit-code correspondence to a particular class of travelling 
wave orbits, namely the fronts. 

Definition 5.1 A front with velocity v is an orbit in A given by 

xl = (t){s -vt- a) Vs, t ez, 

where a & R is fixed and, the front shape (p : M — > R, is a right continuous 
function which obey the conditions 

(f){x) < c if X < 

0(2;) ^ c if X ^ 0. ^ ^ 



The front shape has the following spatio-temporal behavior: 
lim (p(x) = and lim (p(x) = 1. 

In this way, the fronts are actually the travelling interfaces as described in the 
introduction. Moreover, for any front shape, the front changes by varying cr; but 
if V = I with p and q co-prime integers, there are only q different fronts that 
cannot be deduced from one another by space translations. On the other hand, 
when V is irrational, the family of such orbits becomes uncountable. (Both these 
claims are deduced from the proof of Theorem below.) 



The existence of fronts is stated in 

Theorem 5.2 Given (a, e) there is a countable nowhere dense set Ga^e C (0, 1] 
such that, for any c G (0, 1] \ Ga,e, there exists a unique front shape. The cor- 
responding front velocity is a continuous function of the parameters with range 

[-1,1]- 

In other words, for any velocity v, there is a parameter set for which the only 
fronts that exist are those of velocity v. 

Furthermore, for c G Ga,e, no front exists. But, the front velocity can be 
extended to a continuous function for all the values of the parameters. 



Referring to a similar study of a circle's map |15], we point out that, even 
when no front exists, numerical simulations show convergence towards a "ghost" 
front. Actually, by enlarging the linear stability result, this comment is proved 



see Proposition |6]^ below). A ghost front is a front-like sequence in Ai^, that is 
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not an orbit of the CML and, for which the (spatial) shape obeys, instead of 
the conditions 

J 0(a;) c if x <0 
\ (j){x) > c if X ^0. 

Ghosts fronts arise in this model because the local map is discontinuous. More- 
over, given the parameters, one can avoid this ghost orbit by changing the value 
of the local map / at c. 



Proof of Theorem 5. 



It follows from Definition pTT] that if a front with velocity v exists, the code 



associated with this front is given by 

9\ = H{s~vt-a) ys,tEZ, 

where H is the Heaviside function^ Hence by the corresponding front shape 
has the following expression: 

+00 k 

= H H L,kH{x -n + vk) Vx e M. (8) 

k=l n=—k 

Such a function is increasing (strictly if v is irrational), right continuous and its 
discontinuity points are of the form x = n — vk (fc ^ 1, n G Z and \n\ ^ k). Now 
one has to prove that, given (a, e, c), there exists a unique front velocity such 
that (H) satisfy the conditions (|^), i.e. there exists a unique v such that the code 
H{s — vt — a) is admissible. For this let 

+00 [vk] 

r]{a,e,v) = 0(0) = X! ^".fc' 

k=l n=~k 

where the continuous dependence on a and e is included in the ln,k which are 
uniformly summable. It is immediate that 77(0, e, v) is a strictly increasing (e 7^ 0), 
right continuous function of v. Moreover, it is continuous on the irrationals andQ 

r]{a,e,-) -ri{a,e,- ) > 0. 
q q 

The conditions (|^ impose that v be given by 

min {t; G R : iji^a, e, f ) > c} . (9) 
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H{x) 



if X <0 

1 if x^O 



9 „(£-) = lini g(h) 
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Actually, if v{a, e, c) is irrational then 



(j){x) < 0(0) = c if X < 0. 



If v{a, 



e,c) is rational we have 



(j){x) = 0(0 ) = ?7(a,e. 



) if ^ a; < 



Q 



hence either 0(O~) < c and the front exists, or 0(O~) = c and the first condition 
in d^) is not satisfied. In the latter case, given (a, e) and v = ^, there is a unique 
c that realizes <p{0~) = c. The countability of the values of c for which there is 
no front then follows. 

In this way, we can conclude that for a, e and | fixed, there exists an interval 
of the parameter c given by 

p~ p 
V{a,e,- ) < c ^ 77(a,e, -), 
q q 

for which the front shape uniquely exists and the front velocity is v{a,e,c) = ^. 
Moreover, we have 



The continuity of v{a, e, c) is ensured by the following arguments. Given 6 > 0, 
(ID implies that 

?7(a, e, -^(0, e,c) — 6) < c and rji^a, e, v{a, e, c) + 6) > c. 
Then for (a, e, c) in a small neighborhood of (a, e, c), one has 

rji^d, e, v{a, e,c) — 6) < c and rj^d, e, 'u(a, e, c) + 5) > c, 
and again from (^) 

v{a, e,c) — 6 < v{d, e, c) and v{a, e,c) + 6 ^ v{a, i, c), 

hence v{a, e, c) is continuous. The range of v{a, e, c) follows from the continuity 
and the values for e 7^ 0: 



G. 




P 



e (-1,1] HQ 



V{a,e, (-1) ) = 0, T]{a,e, 



1) > and r){a, e, 



1) = 1. 



□ 
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Remark 5.3 Notice that, from ^ and the properties of the function rj, v{a, e, c) 
is an increasing function of c with range [—1, 1] when c varies from to 1 (e ^ 0). 
Further, it has the structure of a Devil's staircase. 

This proof gives also a practical method for computing numerically the front 
velocity by inverse plotting r] versus v (see Figure 1). This picture reveals that 
v{a,l,c) < 1 for some values of c. Notice also that for c ^ ^, the velocity is 
non negative. Moreover, the local map has to be sufficiently non-symmetric, 
i.e. c must be sufficiently different from |, to have travelling fronts. These two 
comments can be stated from the inequality 



6 The linear stability of fronts with rational ve- 
locity 

In this section, we improve the linear stability result of the orbits in A for the 
fronts and, also for the ghost fronts. 

We now consider the configurations x E A4 for which the code is given by the 
Heaviside function, namely the kinks. Among the kinks, we shall use the front's 
configurationsj^: 

+00 k 

= Y. ln,kH{s -n + vk-a) Vs G Z, 

fc=l n=~k 

where v is given by (|^). Further, let the interval 

and [. ] be the ceiling function[^ The linear stability of fronts is claimed in the 
following statement. 

Proposition 6.1 For c G Ip, let S = min |c — 77(|^), 77(2) — c| and x a kink 
initial condition such that 

eV.:\{R^)^^^^^- Xs+\a^\ ^?ia'^ Vfc ^ 0, \J - 2k - I s ^2k, (10) 

then i?^ and have the same code for all times and 

hm =0. 



""^"From now on, the parameters' dependence is removed unless an ambiguity results and, we 
denote by x or by x'^ the initial condition of the orbit {a;*}tgN- 
i.e. [x] e Z and x ^ [a;] < x + \ 
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Remark 6.2 In practice, according to the present phase space, the condition ( |7^ j 
has only to hold for < A; < where to is such that (M + l)a*" < 6. 



Proof: Since x is a kink that satisfies (p^O]), we have according to the nearest 
neighbors couphng of the CML: 



Vt ^ 



^s+M < c if s «c 
^ c if s ^ t. 



Using these inequahties, the condition ([T0| ) and the definition of S, one shows by 
induction, that 

Vt ^ - < 5a"^' V/c^O, V-2A;- 1 -t < s s; 2A; + t. 

The latter induces that both the orbits have the same code for all times. □ 

When it exists, let us define the ghost front of velocity | by {Rpt+a}t&z- For 
such orbits, the linear stability's statement is slightly different than for the fronts: 



Proposition 6.3 For c = //(^ ), let 6 = 0(— ^) — '^'^'^ ^ ^ kink initial 

condition such that 

3a eR: < x^+p^] - iRa)s+ia^ =^ ^^"'^ VA; ^ 0, V - 2A; - 1 ^ s 2fc, 

then 

lim \\Rpt+a ~ x^W = 0. 
The proof is similar to the preceding one noticing that for the ghost front 

Rp(t+l)+a = + ^RPt+a, 

where 9^ = H{s - [|t + a]) and that Vt ^ 



a^s+H > {REt+a)s+\a] Vs. 



7 The interfaces and their velocity 

In this section, we study the dynamics and the properties of the code of the orbits 
for which the state is a kink for all (positive) times. 
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Definition 7.1 An interface is a positive orbit {x*}(gN such that 

^, -1 T/ /\ f X* < c if s ^ J(x^) — 1 

^ ^ [ X* ^ c if s ^ J{x^), 

where the sequence {J{x*)}^^f,J is called the temporal code. The velocity of an 
interface is the limit 

J(x*) 

V = lim , 

provided it exists. 

Notice that the front's temporal code is given by 

\vt + a'] \/te-L, (11) 



namely it is uniform 



To compare the kinks, the following partial order is considered 

Definition 7.2 Let x,y E Ai. We say that x ^ y iff Xs ^ ys \/s E'L. 

As a direct consequence of our CML, using this definition, we have 

Proposition 7.3 (i) If x -< y then F^x -< F^y. 
(a) If X -< y are two kinks, then J{x) ^ J{y)- 

Moreover, we can ensure a kink to be the initial condition of an interface in the 
following ways 

Proposition 7.4 Let x be a kink initial condition. If x^ Xs+i Vs G Z or, if 
e ^ |, then {x*}tgN is an interface. 

We will also use the asymptotic behaviour in space of the interfaces' states: 
Proposition 7.5 For an interface: 

^ / \xi\ ^ a^M if sscJ(x°)-t-l 
^ \\l-xl\ ^ a'M if s ^ J(a;0) + 1. 

This result simply follows from the linearity of the CML and the inequality 

\J{x')-J{x^)\^t, (12) 
which is a consequence of the nearest neighbors' coupling. 

Combining the previous results, one can give some bounds for the temporal 
code of an interface 
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Proposition 7.6 If c&l% and {x*}(j.o is an interface, then 

3^1,^2 GM : Vt ^ r-t + (Til J(x*) ^ r-t + (Tal . 

q q 

P 

Proof: We here prove the left inequahty, the right one follows similarly. If - = —1, 

_ Q 
then the statement holds by the inequality ([1^). Now, let to be given by Remark 
p 

For - 7^ —1 we have 

Let ti be such that a*^M < {Ro)-2to-i ~ J{x^) — ti — 2to — 1- Define 

a;*/ if s < (Ti - 2to - 1 

y,= '{ {Rar)s if ^1 - 2to - 1 < s ^ ai + 2to 

max ^x^f^^ti),xl^ I if s > ai + 2to. 

Let us check that x*^ -< y. For 5"i — 2to — 1 ^ s cti + 2to, ?/s > (-Ro)_2to-i ^ 
using the monotony of R^^ and the previous proposition. 

Hence, by Proposition |775| , we obtain J{x^~^^^) ^ J{y^) Vt ^ 0. But y is a 
kink that satisfies the condition ([T0| ) in the framework of Remark |6.2| , therefore 
the left inequality is stated. □ 

Finally, the statement on the velocity of any interface, valid for all the values 
of the parameters, is given by 

Theorem 7.7 If {x^} ten is an interface then, its velocity exists and 

hm — — = v(a,e,c). 

Proof: For any v given by (D, let — — ^ v, Ci G I^i and Ca G 1% 

qi q2 91 92 

respectively. If c G Ip, then the statement clearly follows from the previous 

proposition. Moreover, by the monotony of the local map, we have 

\/x e M F^^ciX -< F^^cX -< F^^dX, 
where the dependence on c has been added to the CML. It clearly follows that 

Pi . rAKc^^y AKc^) , P2 

— hmmt ■ < limsup ■ < — . 

qi t^+oo t t^+oo t q2 

The convergence is then ensured by choosing — and — arbitrarily close to v. □ 
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To conclude this section, we mention that all these results and, in particular 
Theorem |7.7| , can be extended to the orbits for which the initial condition satisfies 

lim H(xs) - H(s) = 0. 

s— >±oo 



For such orbits, one can prove Theorem |7.7| for the following codes: 

Jinf{x^) = min{s : x* ^ c} and Jsup{x^) = max{s : < c}. 
Actually, we have 

Jinfi-^ ) 1. Jsupi-^ ) — 

lim — — — - = lim — — — - = V. 

8 The nonlinear stability of fronts with rational 
velocity 

An extension of the linear stability result for the fronts can be achieved, also 
using the CML contracting property which appears in (P), by proving that for 
all the kink initial conditions, the corresponding interface's code is identical to a 
front's code, at least after a transient. 



8.1 The co-existence of fronts and non-uniform fronts 

We shall see below that one can only claim the nonlinear stability of fronts in 
some regions of the parameters. We now justify such restrictions on the latter. 

Let, a non-uniform front, be an interface in A for which the temporal code 
cannot be written with the relation (pJ]). Using the admissibility condition, one 
can show the existence of such orbits given a non-uniform code. In this way, we 
prove in Appendix the existence of the non-uniform fronts with velocity | and 
respectively, when e is close to 1. Since these orbits are interfaces, they co-exist 
with the fronts of the same velocity by Theorem \7.1\ . Moreover, it is possible 
to state the linear stability of these orbits similarly to Proposition Hence, 
when the non-uniform fronts exist, the fronts do not attract all the kink initial 
conditions. 



8.2 The nonUnear stabiUty for weak coupUng 

Here, we state the nonlinear stability's result for the fronts with non negative 
rational velocity, in the neighborhood of e = 0. However, the following assertion 
can similarly be extended to the fronts with negative velocity. 
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Let ^ 

k=l 



1 ' q' q 



and, for 6 G (0, 1) let the interval 



we have: 



Theorem 8.1 Given | G Q n [0, 1], a e (0, 1) and 9 G (0, 1), there exists eo > 
such that for any kink initial condition x: 

VeG(0,eo) Vc G li 3a G R : lim - x*]] = 0. 

q t — ^ + 00 



Remark 8.2 For the velocity ^ = 0, the statement still hold for 6* = 0. 

Moreover for the velocity | = 1, using similar techniques , we have proved that 
the theorem holds with ^ = and eo = 1 for all the interfaces. 



8.3 Proof of the nonlinear stability 



In the proof of Theorem we assume e < g then, by Proposition |7.4| , all the 
orbits {x*}f6N under consideration are interfaces. 

In a first place, it is shown that the interfaces propagate forward provided 
some conditions on the parameters are satisfied. 

Proposition 8.3 It exists ei > such that for any interface {x*}tm o-nd any 
\, we have 

Ve G (0,ei) 3to : Vt ^ to ^ J{x'). 

Proof: According to the present phase space, all the orbits of the CML are 
bounded in the following way: 

V5 > 3to : yt^to -S <xl<l + 6 Vs G Z. 

In particular, for an interface: 



<a{l- -)€+- + — \/t^ to. 



Hence, if e is sufficiently small, the statement holds for all c ^ |. □ 
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Furthermore, we introduce a convenient configuration. Then, after computing 
the code of the corresponding orbit, we show a dynamical relation between the 
code of any interface and the code of this particular orbit. 



Let 

+00 k 

\/j e Z (Sj)^ = E E ln,kH{s -n~j) Vs G Z. 

k=l n=—k 

Notice that if c G {r]{Q^),r]{Q)], Sj is the front of velocity 0. When c ^ (??(0~), r7(0)], 
this configuration is a convenient tool. Actually, the temporal code for the orbit 
{'S'ojtgN is shown to be given by: 

Proposition 8.4 Given | > 0, a G (0, 1) and 9 G [0, 1), there exists eg > such 
that 

Ve G (0, eo) Vc G l| J{Sl) = \{t + 1)^] Vt ^ 0. 

q 

Moreover 

inf |(5'*),-c| >0. 
The proof is given in Appendix 0. 

We now state the main property of an interface's code. 

Proposition 8.5 For an interface we have, under the same conditions of the 
previous proposition 

3no Vn ^ no J(x*+") J{S'o'^) + J(x") Vt ^ 1. (13) 

Proof: Using the relation for t = n — 1, we obtain the relation between the 
interfaces' states and the configuration Sj: 



By induction, the latter leads to 

n-l 



fc=m+l 



Then the monotony of Sj, the positivity of L and Proposition ^.3| induce the 
following inequality for m large enough and n > m 



where 



n 
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Consequently 

J(x*+") J{Si) Vt > 0. 



Now, given 5 = inf \ {Sq)s — c\ according to the previous Proposition, let no > 
m be such that 

It follows that 

\\Sn - Sj(^n-i)\\ < 6 Vn ^ no + 1, 
and hence, by Proposition |4.1| 

= J(S'^) + J(x"-i) Vt>0,n^no + l. 

□ 

Finally, we state the result which ensure a sequence to be uniform. 

Lemma 8.6 Let {j„}„gN be an integer sequence which satisfies the conditions 
Vn ^ no.' 

(i) Jn+k ^ \'^k] +jn \/k ^ 0, 

(a) in ^ + a] for a fixed cr G R, 
then 

P 

37 G M and ni > no : Vn ^ ni j„ = T-n + 7] . 

q 

The proof is given in Appendix 0. 

Collecting the previous results, we can conclude that, under the conditions of 



Proposition ^]3| and p.4| , the temporal code of any interface satisfies the condition 
(i) of Lemma |8.6| (Notice that Proposition |8.3| only serves in the proof of Propo- 
sition |8.5|). Moreover, this code also satisfies the condition (ii) by Proposition 



|7.6| . Hence, after a transient, all the interfaces have a front's code. By (y), they 
consequently converge to a front. 



To conclude the proof of Theorem in the framework of Remark |8.2| , we let 
c e Iq. Hence ^(5*0) = Vt. If c > |, the statement holds using Proposition ^ 
and the relation (|T3|). For c < |, one can prove the following inequalities: 

3to : Vt ^ to J{x'^^) ^ J{x'), 

similarly to Proposition ^.3| and, 

3ni : Vn^ni J(a;''+^) ^ 7(50) + J(x"), 



similarly to the relation (TTBI). 
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9 Concluding Remarks 



In this article, we have considered a simple space-time discrete model for the dy- 
namics of various nonlinear extended systems out of equilibrium. The (piece-wise) 
linearity of this bistable CML allowed the construction of a bijection between the 
set of global orbits and the set of admissible codes. When they do not reach the 
discontinuity, these orbits are linearly stable. For e < |, the CML is injective on 
Ai, then A can be identified with the limit set 

+00 

t=0 

which attracts all the orbits. These comments justify the study of the global orbits 
and in particular, the study of fronts which occur widely in extended systems. 

The existence of fronts, with a parametric dependence for their velocity, has 
been proved using the spatio-temporal coding. The velocity was shown to be 
increasing with c. We have in addition checked numerically that v is also an 
increasing function of e on (0, |) for any a and on (0, 1) for any a G (0, |) or for 
any a if c is such that v{a, 1, c) = 1. Moreover, one can find some values of a and 
c for which the front velocity does not (always) increase with e (see Figure 2). 
The spatio-temporal coding also serves to show the existence of other patterns 
such as non-uniformly propagating interfaces (Appendix 0). When these exist, 
they always co-exist with the fronts of the same velocity. 

Furthermore, we have consider the more general dynamics of interfaces. Using 
the temporal code, we have shown that all these orbits have the front's velocity 
uniquely determined by the CML parameters. 

The stability of fronts was also proved, firstly with respect to initial conditions 
close to a front state in their " center" , and secondly with respect to any kinks for 
the fronts with non negative rational velocity, assuming some restrictions on the 
parameters. Actually, the latter allows us to avoid the existence of non-uniform 
fronts which would attract certain interfaces. 

Finally, notice that all the results stated for the orbits that never reach the 
discontinuity can be extended to some CML's with a C°° bistable local map. 
Actually, one can modify the local map into a C°° one, in the intervals where the 
orbits in A never go, without changing these orbits. In other terms, all our results 
stated in open sets can be extended to differentiable maps, in particular, the 
existence and the linear stability of fronts and non-uniform fronts with rational 
velocity. 



17 



These last results show the robustness of fronts in these models. This empha- 
sizes the reahstic nature of these models. 
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A Properties of the coefficients In^k 

In first place, if /c > 1 tlie following expression holds 



ln,k = (1 - a)a'-\-ir E 

p=\n\ 



2 J \pj \p + nj 
where stands for the binomial coefficients. In particular, one has 



lim/„,fc = (1 - a)a'' ^Sn,Q, 



and 

l—a/a\kf k \ -r k+n ^ rj, 

if ^ ^ z 



i'n,k 



Moreover, one can show that these coefficients have the following generating 
function 

— + a(l - e)x + —x^ = V /n-fc,fca^ e M, k ^ 1. 

a \2 2 y 

Prom this function we deduce the relation 

E ^n-m,tlm,k = Vt, A; ^ 1. (^1) 

Finally, we prove the following behaviour when e is small 

Proposition A.l Vg > 1 and q ^ 

(i) lim^ > 0. 

2 +00 k 

(^^^ ^I^T^ ^ ^ ^^^^ = 
A;=5+l n=p+l 

Proof: Assume q ^ I and p ^ 0. (i), From the explicit expression of the Z^^^, we 
have 

^.^ ^ (1 - a)a^-^ /g\ 
6^0 eP 2P \pj' 

Moreover, (ii) is ensured by the relation: 

□ 
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B The existence of non-uniform fronts 

B.l The non-uniform front with velocity | 

Consider the following sequence: 

and 

+00 k 

^1^12 12 kkH{s -n- jt-k) ys,tez 

k=l n=—k 

the corresponding orbit in A, namely a non-uniform front with velocity ^. Actu- 
ally, jt+4 — jt + '2 but it does not correspond to the temporal code of the front of 
velocity |. The existence of this particular orbit is stated in 

Proposition B.l 3e > such that Ve e (e, 1], there exists an interval I so that 

Vcei Vtez = 

Proof: z\ is monotone in s. Then one has to prove that, in a neighborhood of 
e=l, 

min 2;* > max 2* , 
t t J* ^' 

in order for the statement to hold for I = ( max^;*- 1 , min^;! . 

\ t ^* t ^\ 

Let = jt - jt-k- We have 

A*+^ = A* and A*+4 = A*+2, 

and 

minAl = A? and maxAl = A9. 

^ ft (V ^ /v rv 

Using these properties, we show that 

+00 

min 4 - max4_i = 5I(-l)''^fc,2fc- 

k=l 

Prom the expression of the /fc,2fc for e = 1, it follows that 

lim (min 2;! — max 2;! 1 ) > 0. 

e-*l- t 3* t ^* ^' 

□ 

Notice that (when it exists) this non-uniform front lies between the fronts of code 
L^Jand LIJ. 
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B.2 The oscillating front 

Here, we consider a non-uniform front of velocity 0. In particular, the one for 
which the temporal code writes: 

2 

Naturally, the expression for this orbit can be written similarly to the previous one 
but now, with this code's relation. Let {Z^}t& denotes the orbit. Its existence 
is claimed in 

Proposition B.2 3e > such that Ve G (e, 1], there exists an interval I so that 

Vcel VtGZ F,(Z*) = Z*+^ 

Similarly to the previous proof, the existence of this orbit is ensured by the 
relation 

+ 00 

niinZ*^-maxZ*^.i = 5:(-l)%,,. 

k=l 

Furthermore, using these technique, we can prove the existence of another 
non- uniform fronts. All these orbits exist only in a neighborhood of e = 1, they 
disappear for e close to 0. 



C Proof of Proposition 8.4 

Since the fronts' configurations are increasing functions in space, the following 
ordering holds: 

Sq -< Rv. 



Hence, using Proposition |7T3| , we obtain for c G Ip : 

^ °n|(m)i-i " V [£(*+!)] -1 

To prove the statement, we then show by induction that 

in a neighborhood of e = 0. 

Using the relations (H) and (Al) to write 

t k +00 k 

('^o), = E E ln,kH{s-n-J{Si-'))+ E ln,kH{s-n), 

k=l n=—k k=t+l n=—k 



21 



we consider the following decomposition 
where 

Pi = i:j:in,k{H{\^^{t + l)]-\^^it-k+l)]-n)-H{^^k-n] 

k=l n=—k 

t+1 

P2 = Y: ln,t+l{H{\^S^ + l)^-n)-H{^{t + l)-n] 

n=-t-l 

^3 = E j:in,k{Hi\^^it + l)]-n)-Hi^^k-n)) 

k=t+2 n=-k 
+00 

Pi = dy^jkpM- 

k=l 

We have the following inequalities: 
Pi ^ since [2(t + l)] - [2(t-A; + l)] > [£A;J VA;,teN. 
P2 ^ and F2 = ^ [£(*+!)], t+i if t<q-l. 

+00 k 

P3 > i: j:kk{Hi\^^{t+l)]-n)~H{k-n)) 

k=t+2 n=-k 

+00 k 

k=t+2 n=rE(t+l)]+l 

Then one can write, for < t < — 1: 

+00 k 

Pl+P2 + P3 + P4^llE(t+l)],t+l- E E 

" fc=i+2n= [2(4+1)]+! 



and for t ^ g — 1: 



+00 k 

Pi + P2 + P3 + Pa^01p,,- E E 

A;=ij+1 n=p+l 



Hence using Proposition |A.1| , in both cases Pi + P2 + P3 + P4 > if e is sufficiently 
small. □ 



D Proof of Lemma 8.6 



By (n), it exists rii ^ no such that 

p . p 



q q 
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Then, also using (i), one has 

q q Q 

which imphes that 

P 

Now it is straightforward to see that 

P 

in = \n- +7] Vn > rii, 
Q 

where 7 = jn^ - 
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Figure 1: Grey- level plot of the front velocity versus e (Horizontal axis) and 
c (Vertical axis) for a = 0.6 in range (white) to 1 (black). This figure is a 
projection of a two dimensional Devil's staircase. 
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Figure 2: Plot of the front velocity versus e for a = 0.9 and c = 0.8216. 
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